Two-finger selection theory in the Saffman- Taylor problem 
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We find that solvability theory selects a set of stationary solutions of the Saffman- Taylor problem 
with coexistence of two unequal fingers advancing with the same velocity but with different relative 
widths Ai and A2 and different tip positions. For vanishingly small dimensionless surface tension do, 
an infinite discrete set of values of the total filling fraction A = Ai + A2 and of the relative individual 
finger width p — Ai/A are selected out of a two-parameter continuous degeneracy. They scale as 



A — 1/2 ~ d and \p — 1/2 1 ~ d ' . The selected values of A differ from those of the single finger 
case. Explicit approximate expressions for both spectra are given. 

PACS numbers: 47.54. +r, 47.20. Ma, 47.20.Ky, 47.20.Hw 



The Saffman- Taylor problem has played a central role 
in the field of interfacial pattern selection in the last few 
decades [1]. It deals with the morphological instability 
of the interface between two inmiscible fluids confined in 
a quasi twodimcnsional (Hcle-Shaw) cell, when the less 
viscous fluid is displacing the more viscous one in a chan- 
nel geometry [2]. In particular, in their seminal work [3] 
Saffman and Taylor called the attention upon the so- 
called selection problem, namely the fact that a unique 
finger-like steady state solution is observed whereas a 
continuum of solutions is possible if surface tension is 
neglected. Full analytical understanding of the subtle 
role of surface tension acting as the relevant selection 
mechanism was not achieved until much more recently 
[4-8] . The resulting scenario of selection has been shown 
to apply with some genericity to other interfacial pattern 
forming systems, most remarkably in dendritic growth 
[1,9]. On the other hand, despite the relative analytical 
tractability of the problem, the dynamics of competing 
fingers is far from being understood even at a qualitative 
level. Recently, it has been shown that in general surface 
tension may affect the long time dynamics in an essential 
way [10]. In the case of the dynamics of finger arrays, the 
effect of surface tension becomes particularly dramatic, 
showing that the qualitative picture of finger competition 
based solely on the concept of screening of the laplacian 
field or the global instability [11] of a periodic finger array 
is insufficient [12]. 

Existence of multifinger stationary solutions of the zero 
surface tension problem has been known for a long time. 
In Ref. [12] it has been emphasized that multifinger sta- 
tionary solutions are relevant to the issue of the dynam- 
ical role of surface tension. In particular the equal-finger 
fixed point has been pointed out as the relevant saddle- 
point to describe competition dynamics. In connection 
with the phase flow structure around this fixed point, 
the problem of existence of unequal-finger fixed points 
with nonzero surface tension has been posed [12]. Here 
we will extend selection theory to search for such solu- 
tions in the case of two fingers. We will follow the for- 
mulation of Hong and Langer [5], which is based on a 



Fredholm solvability analysis of a non-self adjoint prob- 
lem defined through linearization about the zero surface 
tension solution, together with WKB and steepest de- 
scent techniques. Despite the admitted objections to the 
full quantitative validity of the approximations involved 
in this method [6,8], it has been shown that it leads to 
the correct qualitative picture of selection and the correct 
scaling of solutions [8]. It is therefore suitable, for sim- 
plicity of calculus and presentation, for a first exploration 
of novel situations such as the present one. 

Our starting point is the dynamical equation for the 
conformal mapping f(w, t) which maps the interior of the 
unit circle in the complex plane w into the viscous fluid 
region, with the unit circle w = e 1 ^ being mapped into 
the interface. Without loss of generality, we will assume 
a channel width W = 2tt in the y direction (with peri- 
odic boundary conditions) and a velocity Uqo = 1 of the 
fluid at infinity. We define the velocity of the stationary 
solutions of the interface as U = j where A is the total 
filling fraction of the channel by the invading fluid. The 
cartesian coordinates in the frame moving with velocity 
U are given by z — x + iy — f(w, t) — Ut. The mapping 
f(w, t) contains a logarithmic singularity which is due to 
the fact that we are mapping an unbounded domain (the 
scmiinfinite strip) into the unit circle, in such a way that 
f(w, t) + \nw is always an analytic funtion in the interior 
of the unit circle. The exact dynamical equation for the 
mapping can be written in the form 



i?e(zd /(rM)dJ*(<M)) = 1 - Udod+HtiK] 



(1) 



which can be easily derived for instance from Ref. [2]. 

Here do is a dimensionless surface tension defined as 
,2 

do = ™ u , where a is the surface tension, b is the gap 
and p, is the viscosity. The curvature k can be expressed 
in terms of x((j>) and y(4>) as 



djxd^y - djyd^x 



(2) 



H is a linear integral operator (Hilbert transform) which 



1 



acts on a real 27r-periodic function g(<j>) according to the 
definition 



1 f 27T 1 
H 4>[9\ = T^P J o 9(s)cotg-(4> - s)ds. 



(3) 



It follows from Eq.(3) that the function A(w) defined by 
A(e 1 ^) = g((j>) + iH^[g] is analytic in the interior of the 
unit circle, and Im(A(0)) = 0. 

In the steady state we will have dtf*((f>,t) = U, and 
Eq.(l) will read 



1 - Ud 



d<p 



(4) 



We search for solutions of the generic type described in 
Fig.l. The total filling fraction A is splitted into two 
contributions Ai + A 2 = A, and we define as a new selec- 
tion parameter the relative finger width p = Xi/X. For 
simplicity we will consider fingers which are axisymmet- 
ric and for convenience we will fix the tip positions at 
(f> = 7r/2,37r/2, for all A and p. The filling fraction A 
ranges from to 1. We take A2 < Ai so that p ranges 
from 1/2 to 1. In these conditions the two fingers corre- 
spond to the intervals 0i = ^(1 — 2p) to 2 = §(1 + 2p) 
and 02 to 3 = 2n + 0i respectively. 

After integration over of Eq. (4) we obtain 



Uy(<f>) 



UdoH^n] + c(0) 



(5) 



where c(0) is a piecewise constant function. The values 
it takes at the intervals (0i, 2 ) and (02, 03 ) differ by the 
finite amount tt(U — 1) which accounts for the disconti- 
nuity of y and the finite flux at the points 0i , 02 • 

After Hilbert transform of Eq.(5) and using that 
H l\9] = -5(0), we obtain 



- d O K(0) + x(0) = X O (0) 
y(0) + = const + H^[x\. 



(6) 
(7) 



where Eq.(7) is just an expression of analyticity of 
f(w,t) + Inw. The function x o (0) is found explicitely 
as Xo(0) = Hj,[g] with <?(0) = (A — 1)0 + c(0), and by 
construction it corresponds to the solution of the zero 
surface tension case. In our case it reads 



xq(0) = (1 — A) ln(2| sin0 — cosp7r|). 



(8) 



Completed with yo(0) = — M0 + c (0))> this gives a two 
parameter class of exact solutions of the type of Fig.l, 
for do = 0. Both A and p can be varied continuously 
within their natural range. The difference between 
the x-coordinate of the two tips is given by 



(1-A) In 



1 — COSp7T 
1 + COSP7T 



(9) 



These solutions are precisely those studied in Ref. [12]. 

The present formulation has some interesting advan- 
tages over the traditional approach of McLean and 



Saffman [4] for instance in that the zeroth order solu- 
tion is obtained naturally as an explicit outcome of the 
method and that it is more amenable to generalization, 
for instance to a larger number of fingers [13]. 

We now proceed by assuming x(0) = Xo((f>) + <ioXi(0) 
and linearizing on xi(0) but keeping all singular terms 
necessary for selection. (Nonlinear effects are expected 
to introduce only a slight quantitative correction to the 
final spectrum of selection [6,8].) Using the relation 
2/1 (0) = Hc&xi] we get 



, d 2 x\ , , ,. 



d 2 H, p [x l } 



+ r(0)x 1= M0) (10) 



where r(0) and p(0) are given by 



r(0) 



A 2 «[( 5 (0)) 2 + ^cos^], 



9 4 (0) 



P(0) = 



1 cos< 



/3<Z(0) 



(11) 
(12) 



with q{4>) = sin0 — cosp7r and (3 = y^j. Explicit knowl- 
edge of /Lt(0) is not necessary for the solvability analysis. 
First order derivatives are subdominant as do — > and 
have been omitted [5]. 

The linear operator on the lhs of Eq.(10) can be seen 
as a 2 x 2 matrix operator acting on a vector of two 
components x+(0) and £7/(0) which are defined respec- 
tively on the intervals (0i, 02) and (0i, 03). Inserting an 
ansatz of WKB form with a point of stationary phase 
of the solution in the upper (or lower) complex plane [5] 
one can show, using steepest descent techniques, that the 
off-diagonal terms of Eq.(10) lead to exponentially small 
contributions. As a consequence, to leading order the 
problem is decoupled into two separate problems defined 
in two disjoined intervals. Similarly, neglecting exponen- 
tially small terms, the integral part of the diagonal terms 
takes a purely differential form in the complex plane [5]. 
The change of variables n — — ^ sin ™ c ^ s p7r maps sepa- 
rately each of the two disjoint intervals above into the 
whole real axis 77 e (—00,00). Therefore, to leading or- 
der we endup with two (complex) differential equations 
of the form 



d °~d^~ + Q± ^ x ^ = R± ^- 



(13) 



which are mutually independent but linked through the 
dependence on the parameters A and p. More details 
of this derivation will be presented elsewhere [13]. We 
define two solvability functions as 



A±(\,p;d ) 



f 

J — ( 



(14) 



where x (77) are eigenfunctions of the null space of the 
adjoint operators of the respective homogeneous equa- 
tions [13]. To enforce solvability we now have to impose 



2 



the simultaneous vanishing of the two solvability func- 
tions A±(X,p; do) = 0. These two conditions will fix the 
discrete spectra of possible values of both A and p. 

Within the WKB approximation, the two solvability 
functions take the form 



/CO 1 
-oo 



*±(tj) 



d?7 



(15) 



where 



IT 



J" 

Jo 



1 + /3 2 t/ 2 



C0Sp7T 



\A + rj 12 (3 2 sin 2 pn 



dr]'. 



(16) 



In order to estimate the solvability functions in the steep- 
est descent approximation, only the form of ^±{rf) is re- 
quired. The singularity structure of *&±(r)) is such that 
the cases p = 1/2 and p ^ 1/2 must be treated sepa- 
rately. The first case (two identical fingers) degenerates 
into the usual single finger problem. For p > i, a more 
complicated singularity structure is revealed. In the up- 
per half complex plane of 77, we find that ^±isl has a 
new branch point at 77 = i/(3 sinp7r, in addition to the sin- 
gularities that were present in the single finger problem, 
namely, a branch point at 77 = i and a pole at i/ (3. On 

the other hand, ^ has the branch point at 77 = i and 
the new one at 77 = i/(3ampTr 7 whereas the pole at i/(3 
is suppressed. Since l//3sinp7r > 1/(3, we obtain that 
(3 > 1 is a necessary condition for the first solvability 
function A + (A,p;d ) to oscillate, and therefore generate 
zeroes. We thus recover the condition A > 1/2 of the 
single finger case, but now for the total filling fraction. 
The equivalent condition for A_(A,;j;do) is /3sinp7r > 1 
so that the new singularity at 77 = i//3s'mpir stands be- 
low 77 = i. This condition also implies that in the contour 
integration for A + (A,p; do) we will always pick up a con- 
tribution from this new singularity. 

By deforming the contour integral as indicated in Fig. 2, 
and following Ref. [5] in identifying the crossover from 
oscillating to non oscillating behaviour of the solvability 
functions, we obtain the scaling of both A and p with 



do to be (A - i) ~ d 2 / 3 and \p — i| ~ d^ 3 . According 
to Eq.(9), the resulting scaling for the tip difference is 
A x ~ d Q ' . 

An explicit (approximate) discrete spectra of selected 
values of A and p for small do will be given by the condi- 



tion cos( 



^±(i+0)-^±(i-0) - 



0. 



From the condition A_(A,p; do) 



we thus obtain 



1 (1-A) 5 



\fd~o~ ^ 
with to = 1, 2, ... and where 



to 



1 r 

-27 C ° tg Wo 



u*H(u;0,p) 



(u 3 - u)(u! - 11)2 



(17) 



r dti (18) 



with the regular part of the integrand H(u;[3,p) = 
(2 — u)*(u4 — 7i) _1 (7i2 — u)~2 and 7i 1-2 = 1 T l//3sinp7r, 
"3,4 = 1 T !//?• 

Finally, from condition A + (A,;j;do) = 0, expressing 
Eq.(17) to leading order and using properties of hyper- 
geometric functions, the two selection conditions can be 
combined to read [13] 



V3o= (2A - 1)3 = n 
S(a) = I(m-i) 



(19) 
(20) 



with n = 1, 2, ... and where 



of s 3V2tt 5 19 

s ( a ) = 2F1 I -, -; -; 1 



4' 2' 4' 



a (l-a)3. (21) 



2F1 is a hypergeometric function [14] and a = ^- 2A _ 1 
is of order (do) and ranges from to 1. 

Eq.(19) determines a set of discrete values of A. No- 
tice that these are given independently of p but the set 
of values are inserted between those of the single finger 
case (p = 1/2), which in the same approximation arc 
given by — =(2A — l)i = n — \ in place of Eq.(19). On 

the other hand, the lhs of Eq.(20) is a monotonically de- 
creasing function of a which varies continuously from 1 /4 
(at a = 0) to (at a = 1) [14]. Solving Eq.(20) for a pro- 
duces solutions with p / 1/2. These will exist whenever 
— (m —5) < I- For a given n, the solutions are labeled 
by 777 = 1 , 2, ... up to the integer part of (n+l)/4. There- 
fore, the first solution with p ^ 1/2 will appear at n = 3 

and gives \p — || ~ 0.3886 dj 1 / 3 + ... For fixed to, p is an 
increasing function of n (like A), but for fixed n, p has 
its maximum value at to = 1 and then decreases with m 
(See Fig. 3). The spectra here derived must be taken with 
some caution, since they are only approximate. An exact 
calculation to lowest order in d should include nonlinear 
effects and a proper treatment of the turning points in 
the WKB analysis, but the corrections are expected to be 
quantitatively small [15]. More details will be presented 
elsewhere [13]. 

Concerning the stability of these solutions, it is rea- 
sonable to presume that, in general, they will be globally 
unstable, such as established numerically for the equal 
finger array (p = 1/2) in Ref. [11]. This implies that 
they would only be directly observable as a transient 
slowing down of the competition dynamics whenever an 
initial condition is prepared close to any of those solu- 
tions. From a dynamical systems point of view, we would 
like to emphasize that the knowledge of the fixed points, 
even if unstable, is always relevant to elucidate the topo- 
logical features of the phase space flow, and therefore to 
gain insight and qualitative understanding of the dynam- 
ics. In particular the location of these new fixed points 
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We conclude by remarking that, although the present 
solvability analysis is not a rigorous proof of existence of 
solutions it reveals by itself a quite unexpected richness of 
the problem. It would be interesting to search for these 
solutions numerically or by other more rigorous means 
[6,8]. The sole existence of the predicted solutions and 
its presumable generalization to a larger number of fin- 
gers has important consequences on the physical picture 
of finger competition, which turns out to be much more 
complex than common arguments of laplacian screen- 
ing seem to suggest. The common picture, according 
to which fingers slightly ahead escape from their neigh- 
bors, is not necessarily valid in general because of the 
existence of growth modes with unequal noncompetiting 
fingers. For vanishingly small surface tension, however, 
these modes collapse and only the equal-finger multifin- 
gcr mode (p = 1/2) survives as a stationary state. Fi- 
nally, given the genericity of the solvability mechanism of 
selection, this opens the possibility of finding similar so- 
lutions in related problems such as needle crystal growth 
or viscous fingering in circular geometry. 
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Fig. 1. Typical configuration of a two- finger stationary 
solution. 

Fig. 2. Deformation of the steepest descent contour of 
integration in the complex r\ plane with (3 > 1 and p > | 
(a) for A + ; (b) for A_. 

Fig. 3. Spectrum of p as a function of n for different 
values of m. 
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